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Линейное масштабно-пространственное 
представление изображения с помощью 
вейвлет-преобразования 


Разработан метод линейного масштабно-пространственного представления изображения на основе 
интегро-дифференциального уравнения 1-го порядка со сверткой вейвлет-функций в качестве ядра 
интегрального оператора. 


Введение 


Масштабно-пространственное представление — это множество копий исходного 
изображения с разными масштабами, которые позволяют охарактеризовать объекты 
иерархической структуры [1], т.е. объекты, которые содержат подобъекты. При реше- 
нии прикладных задач используются те копии исходного изображения, которые имеют 
наибольшую семантическую значимость в зависимости от цели обработки. 

Масштабно-пространственное представление изображения широко используется 
при обработке медицинских изображений, изображений аэрофотосъемки, а также при 
контроле качества изделий в промышленности, когда решаются задачи определения 
формы объекта, контурной сегментации, определения оптического потока и 
стереозрения. Такое представление позволяет существенно повысить оперативность 
обработки и достоверность распознавания объектов на изображениях. 

Масштабно-пространственное представление изображения может быть линей- 
ным или нелинейным. Линейное масштабно-пространственное представление строится 
на основе уравнения диффузии и ставит в соответствие изображению Гх,у) зависящее 
от масштаба { однопараметрическое семейство изображений Г(х.у,6), которое удов- 
летворяет уравнению 


Ох, у, 1 
9 Ак, у, (0 
01 
2 д? 
с начальным условием Г[(х,у,0) = Цх,у), где А = = + 57’ х, у — пространственные 
х 


координаты. Предполагается, что масштабный параметр Е > 0 монотонно возрастает 
и некоторым образом может быть соотнесен с пространственным масштабом [1]. 
В [1] показано, что Гх,у4) может быть получено с помощью свертки с 


22 
х+у 
41 


гауссовской функцией С(х,у,0) = е . В этом случае построение масштабно- 


1 
41 
пространственного представления сводится к построению последовательности сглажен- 
ных изображений. Производные изображений [(х,у,ё) также удовлетворяют уравне- 
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нию диффузии, а производные гауссовских функций генерируют масштабно-прост- 
ранственное представление, подчеркивающее контуры объектов на разных масш- 
табах. 

Процесс получения масштабно-пространственного представления изображения 
на основе уравнения (1) называется изотропной диффузией. Изотропная диффузия 
обеспечивает помехоустойчивость представления путем сглаживания изображения, 
однако в этом случае размываются контуры, которые часто являются наиболее 
информативной частью объекта распознавания. 

Для решения этой проблемы часть информации о контурах вносят в качестве 
нелинейности в дифференциальное уравнение в частных производных. Тогда внутри 
однородных областей изображение сглаживается сильнее, чем вдоль границ облас- 
тей. Полученное в результате решения задачи Коши для нелинейного дифференциаль- 
ного уравнения в частных производных масштабно-пространственное представление 
изображения называется нелинейным, а подход, который применяется для его 
построения — анизотропным диффузионным [2]. Этот подход позволяет повысить 
качество выделения контуров и улучшить результат сегментации. Уравнение анизо- 
тропной диффузии имеет вид: 


Се [ум] 


где Ф» — дивергенция, ее ояь градиент, У =| —,— | 8(%/) — поро- 
0х ду ох ду 


говая функция для перепадов интенсивности изображения такая, что (0) = 1, 
2(х) > 0, Пт 2(х) = 0. Выбор ®(х) производится по аналогии с выбором порога для |У/ | 
х—>0 


в дифференциальных методах контурной сегментации изображений [2]. Недостатком 
анизотропного диффузионного подхода являются большой объем вычислений по 
сравнению с линейной диффузией и низкая помехоустойчивость из-за использо- 
вания производных по пространственным координатам. Значительная часть практи- 
ческих задач характеризуется высоким уровнем помех. Поэтому задача разработки 
нового подхода к построению масштабно-пространственного представления изобра- 
жения, позволяющего локализовать контуры и обеспечивающего высокую помехо- 
устойчивость, является актуальной. 

В случае иерархической структуры объекта распознавания для получения 
масштабно-пространственного представления изображения вместо пространствен- 
ного дифференцирования в данной работе предлагается применять репагулярное 
вейвлет-преобразование. Это преобразование определено в [3] как свертка строки 


ИГ $2 Хх, х За 


0, 


фиксированное число, зависящее от аЕ(0,1) — параметра преобразования. Применение 
репагулярного вейвлет-преобразования для подчеркивания перепадов интенсивности 
изображения в задаче контурной сегментации изображений позволило достичь более 
высокой помехоустойчивости методов сегментации по сравнению с применением 
дифференцирования [3]. 

Для характеристики помехоустойчивости в работе используется показатель 
качества масштабно-пространственного представления, оценивающий близость резуль- 
татов представления для зашумленного и незашумленного изображений, и пока- 


или столбца изображения с функцией у(х,а)= те: &. = 


Е 
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затель эффективности, учитывающий уменьшение энтропии изображения при сни- 
жении уровня помех. Так как применение репагулярного вейвлет-преобразования при 
построении масштабно-пространственного представления обусловлено необходи- 
мостью сохранения границ объектов при подавлении шума на однородных участках 
изображения, целью работы является повышение эффективности линейного 
масштабно-пространственного представления изображения с помощью вейвлет- 
преобразования. Для достижения поставленной цели решены следующие задачи: 

— разработан метод линейного масштабно-пространственного представления 
изображения на основе интегро-дифференциального уравнения 1-го порядка со 
сверткой вейвлет-функций в качестве ядра интегрального оператора; 

— экспериментально исследованы качество и эффективность предложенного метода. 


Построение линейного масштабно-пространственного 
представления изображения с помощью 
вейвлет-преобразования 


Как уже упоминалось, линейное масштабно-пространственное представление 
строится на основе уравнения диффузии (1). Обозначим т-ю строку изображения 
[(х.у,ё) через и(х,д): 

и(х,!) = (х,ут,ц. 

Тогда (1) для т-й строки изображения имеет вид: 


ди(х,г) _ 0?и(х,г) 
[91 О в 


При моделировании изображения, содержащего перепады интенсивности, 
последние описываются функцией Хевисайда. Поэтому в общем случае функция 
и(х,6), удовлетворяющая (2), имеет разрывы 1-го рода и может быть представлена как 
обобщенная функция, зависящая от параметра 1. Дифференцирование в (2) выполняется 
для обобщенной функции [4]. 

Дифференцирование по пространственной координате можно представить в 
виде свертки строки изображения с производной дельта-функции [4], тогда (2) 
имеет вид: 


(2) 


ЯВ = бк) (6х жи(,0). (3) 
01 

При использовании вейвлет-преобразования для подчеркивания контуров 
изображения задача сегментации сводится к задаче определения точек, в которых 
изменяется регулярность функции значений интенсивности изображения в зависимости 
от пространственных координат. Для обнаружения изменений регулярности функции 
целесообразно использовать специальный вид вейвлет-преобразования, характеризу- 
ющийся не изменением параметра масштаба, а изменением показателя регулярности 
вейвлет-функции. Этот вид вейвлет-преобразования, которое называется репагулярным 
(от латинского герагит — сдерживающая преграда), реализуется как свертка 
изображения с вейвлет-функциями разной регулярности, локализованными в одной 
точке. Известно, что подчеркивание перепадов интенсивности изображения с 


помощью вейвлет-функций разной регулярности у(х,а), где ае(0,1), обеспечивает 
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большую по сравнению с операцией дифференцирования помехоустойчивость [3]. 
Поэтому заменим в (3) 6'(х) на функцию вида у(х,а), где ае(0,1) фиксировано, 
получаем 
ди(х, 1) й _ 
01 


Формула (4) представляет собой обыкновенное интегро-дифференциальное 
уравнение 1-го порядка со сверткой вейвлет-функций в качестве ядра интегрального 
оператора [5]. Поставим задачу Коши: найти решение этого уравнения (обобщенную 
функцию, которая зависит от параметра 1), обращающееся при { = 0 в заданную 
обобщенную функцию Хх). Чтобы раскрыть физический смысл (4), учтем, что 
первообразная функция у(х,а) имеет вид: 


ука) * (уса) *и(х,1). (4) 


—а-1 


+С(а), 


1 
ф(ха) = х 
—а+1 
где С(а) - произвольная функция а [4], поэтому выберем С(а) = 0. В этом случае 


ба) = ба) 
1-Х, 


Тогда (4) можно представить как 
пы обед *[ 2 обод | (5) 


Согласно свойству дифференцирования свертки обобщенных функций, фор- 
мула (5) принимает вид: 


ди(х,!) а т 0?и(х, г) 
к И 6 
5, ф(х, а) [ос а) Эр (6) 
В области преобразования Фурье (6) и к следующему виду: 
оо, = (7) 


Формула (7) представляет собой в области преобразования Фурье уравнение 
диффузии в неоднородной среде. Использование такого уравнения предполагает, что 
сглаживание спектральной плотности изображения можно представить как процесс 
диффузии в неоднородной среде. 

Чтобы найти концентрацию вещества в любой точке среды в любой момент 
времени, недостаточно одного уравнения (7). Необходимо знать еще распределение 
концентрации вещества в начальный момент времени — начальное условие: 


и(@,0) = /(@). Таким образом получаем в области преобразования Фурье задачу 
Коши для уравнения диффузии в неоднородной изотропной среде: 


ди(@, И _ 
01 


обод и, 8 


и(@0,0) = Г (6). 
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Решение задачи Коши (8), описывающей преобразование спектральной плотности 
изображения как процесс диффузии в неоднородной среде, получается путем 
применения преобразования Фурье к решению (4) с начальным условием 


и(х,0) = /(х). (9) 


Для нахождения решения обыкновенного интегро-дифференциального уравне- 
ния |-го порядка (4) выполняется преобразование Фурье обеих частей этого уравне- 
ния. Получаем линейное однородное уравнение относительно и(@,а) с переменными 
коэффициентами: 


ди(@,1) 
01 


Общее решение уравнения (10) имеет вид: 


= (и (а) и(®,0. (10) 


се". 


Выполнив обратное преобразование Фурье по переменной ©, получаем 


% (и(в,а)у т 


Е(х,г) = С, [е ео. 


Фундаментальным решением уравнения (4) называется решение Ё(х,0), удов- 
летворяющее начальному условию Е(х,0) = 6(х). Если обобщенная функция Хх) в (9) 


финитна, то решение задачи Коши для уравнения (4) с начальным условием (9) запи- 
сывается в виде: 


и(х,г) = Е(х,1) * КХ). (11) 


Для определения константы Со воспользуемся начальным условием задачи 
Коши для уравнения (4): 


Е(х,0) =6(х). 
В области преобразования Фурье по переменной х получаем 


Е(@,0) =С, =1. 
Тогда 


2 1 
Е(х,г) = ее ‘е“а@, 


причем полученное решение не ограничивает выбор вейвлета только репагулярной 
вейвлет-функцией. 


Докажем, что обобщенная функция и(х,!) = Е(х,6)* Хх) удовлетворяет уравне- 
нию (10). Выполнив преобразование Фурье по переменной х, имеем 


и( 0,1) = Е(®,1) 1 (©). (12) 
Подставляя (12) в уравнение (10), получаем 


се Г) 


Хе) — ив.) Е, = Г оо 2) бодр Евы -0 
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Докажем, что (12) удовлетворяет начальному условию задачи Коши и(х,0) = Хх). 
В области преобразования Фурье по переменной х это условие принимает вид: 


и(©,0) = Г(@). 
Тогда 
Шт и(@, а) = Ша ©" (0) = (6). 


Реализация линейного масштабно-пространственного 
представления изображения с помощью 
вейвлет-преобразования 


Линейное масштабно-пространственное представление изображения с помощью 
вейвлет-преобразования согласно формуле (11) реализуется в виде свертки каждой 


строки, а затем каждого столбца матрицы значений интенсивности изображения С 


ильтром {е (п)! , коэффициенты которого представляют собой дискретные значения 
р ИЛ 50 р р р 


функции Е(х,). Для вычисления коэффициентов фильтра {е, (п) задается параметр 
масштабно-пространственного представления ае(0,1). Затем определяются значения 
функции и (о, а) для последовательности дискретных значений ®, например, из 
интервала [- 3, 3] с шагом 0,25. Для упрощения вычислений в качестве и (о, а) 


использовалась функция й(@,а)= 21с0$ 2 а+1]®|“ зп @, являющаяся пре- 


— —а 
образованием Фурье (ха) = х $21 х. Результат вычисления интеграла для каждого 


значения частоты потенцируется, а затем к полученной последовательности значений 
применяется обратное дискретное преобразование Фурье. Получаем коэффициенты 
№ 
фильтра {© (п) (рис. 1, 2). 
_ т т т т т т | | т а 
9,0 С(®, 8) 
: : : : : : : 0.9 


1 


0.8 


0.6 


0.4 


0.2 


а) 
Рисунок 1 — Функция, производящая масштабно-пространственное 
представление в пространственной (а) и частотной (6) области: 
1 — гауссовская (#= 1/32), 2 — предложенная (#= 1/4) 
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Рисунок 2 — Масштабно-пространственное представление изображения: 
линейное (а) (1=2, 4, 8, 12, 20), линейное с помощью вейвлет-преобразования 
(6) (Е= 1/24, 1/12, 1/6, 1/4, 1/2) 


Экспериментальные исследования и выВОДЫ 


Экспериментальные исследования предложенного масштабно-пространственного 
представления проведены на изображении размером 256 х 256 пикселей, в центре 
которого — белый квадрат 64 х 64 на черном фоне. На это изображение накладывался 
аддитивный гауссовский шум, для которого отношение сигнал/шум 4 определялось 


по формуле а=й’/с”, где й — величина контраста объекта, о — стандартное 
отклонение шума. 
Показателем качества масштабно-пространственного представления, применяемого 


для улучшения изображения, выбрана величина Р’,„, ‚ определяемая по формуле [6]: 


ро › 


) 


где [”(1,/) - результат масштабно-пространственного представления тестового 
изображения; /“ (1, /) — результат масштабно-пространственного представления зашум- 
ленного изображения; /° (1, /) — зашумленное изображение. 

Для оценки эффективности масштабно-пространственного представления 


изображения предложено использовать коэффициент уменьшения энтропии, так как 
снижение уровня помех приводит к уменьшению энтропии изображения [6]: 


вы 
НР 


уй 
где Н — энтропия зашумленного изображения; Н = -У Р(х,) 05, Р(х,); Ё - количество 
11 
значений интенсивности изображения; Р(х,)- вероятность появления значения 
интенсивности х,; НР- энтропия изображения, полученного в результате масштабно- 
пространственного представления. 
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Получены графики зависимости значения показателей эффективности и качества 
масштабно-пространственного представления изображения от отношения сигнал/шум 
по мощности (рис. 3). 


05 
Ро 


2 5 10 20 50 100 2 5 10 20 50 100 
9 по мощности д по мощности 


Рисунок 3 — Зависимости показателей качества (а) и эффективности (6) 
от отношения сигналЛпум по мощности для масштабно-пространственного 
представления: линейного (1), предложенного (2), нелинейного (3) 


Анализируя полученные результаты, заметим, что предложенное масштабно- 
пространственное представление изображений, пораженных аддитивным гауссовским 
шумом, превышает по эффективности линейное и нелинейное масштабно-пространст- 
венное представление в 1,4 — 2,1 раза при значениях отношения сигнал/шум 100 и 
менее по мощности. 

По значениям показателя качества предложенное масштабно-пространственное 
представление хуже линейного представления до 5 раз при отношении сигнал/шум 
20 и менее по мощности. При значениях отношения сигналлпум 50 и более по 
мощности предложенное представление превышает по показателю качества линейное 
масштабно-пространственное представление до 2,3 раза. Сравнение предложенного 
масштабно-пространственного представления с нелинейным представлением показало, 
что оба представления сходны по показателю качества. 

По быстродействию предложенное масштабно-пространственное представление 
сравнимо с линейным представлением и лучше нелинейного представления в Си раз, 
где п — количество итераций при вычислении нелинейного масштабно-пространст- 
венного представления (в наших экспериментах 30 - 45), С — постоянная, которая 
приблизительно равна 0,1 (рис. 4). 
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Рисунок 4 — Время вычисления нелинейного масштабно-пространственного 
представления изображения 256 х 256 пикселей в зависимости от количества 
итераций (вычисление предложенного масштабно-пространственного представления 
заняло 0,5 — 0,6 с) 
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Таким образом предложенное линейное масштабно-пространственное представле- 
ние целесообразно использовать для улучшения изображений, пораженных аддитивным 
гауссовским шумом при отношении сигнал/шум 100 и менее по мощности в задачах, 
где требуется высокое быстродействие и эффективность. 
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Лнне масштабно-просторове представлення зображення за допомогою вейвлет-перетворення 

Розроблено метод лийного масштабно-просторового представлення зображення на основ! 1нтегро- 
диференшиного р1вняння 1-го порядку 3! згорткою вейвлет-функщй у якост! ядра 1нтегрального 
оператора. 
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